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Abstract 

We discuss the relation between density matrices and the uncertainty 
principle; this allows us to justify and explain a recent statement by 
Man'ko et al. We thereafter use Hardy's uncertainty principle to prove a 
new result for Wigner distributions dominated by a Gaussian and we re- 
late this result to the coarse-graining of phase-space by "quantum blobs" . 

1 Introduction 

Identifying the class of all phase-space functions who are Wigner distributions 
of some mixed quantum state is a formidable and unfinished task. The problem 
is actually the following: assume that a function W(x,p) on phase space defines, 
via the Weyl correspondence, a self-adjoint operator p with unit trace. Then p 
is, a priori, an excellent candidate for being the density operator of some mixed 
state with Wigner distribution W(x,p), provided that in addition this operator 
is positive (that is (pip \tp) > for all square integrable ip). And this is where 
the difficulty comes from: outside a few well-known cases (for instance when p 
is a Gaussian), it is notoriously difficult in general to check the positivity of p 
by simply inspecting the function W(x,p) (the condition W(x,p) > is neither 
necessary nor sufficient, in strong opposition to the classical case). Although 
there are general (and difficult) mathematical theorems giving both necessary 
and sufficient conditions for positivity (the "KLM conditions" , which we shortly 
review in Section [2| , these results are not of great help in practice because 
they involve the simultaneous verification of the positivity of infinitely many 
square matrices of increasing dimension. A supplementary difficulty is actually 
lurking in the shadows: these conditions are sensitive to the value of Planck's 
constant when the latter is used as a variable parameter: a given operator p 
might thus very well be positive for one value of h and negative for another 
(this somewhat unexpected but crucial property is best understood in terms of 
the Narcowich- Wigner spectrum [HE])- This feature is of course completely 
fatal when one wants to use semiclassical or WKB methods. In fact, in a recent 
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very interesting Letter [3] Man'ko et al. have shown that rescaling the position 
and momentum coordinates by a common factor can take a density matrix into 
a non-positive operator while preserving a class of sharp uncertainty relations 
(the Robertson-Schrodinger uncertainty principle, which we recall in Section 
[2]). This allows these authors to conclude that "...the uncertainty principle does 
not determine the quantum state" . Man'ko et al. are of course right; in fact 
Narcowich and O'Connell 4J had already shown in the mid 1980s that fulfilling 
the uncertainty relations is necessary, but not sufficient, to ensure the positivity 
of p (this example is described in next Section). 

The goal of this Letter is threefold. First (Section [2]), we complement and 
explain from a somewhat different (and more critical) perspective the results 
of Man'ko et al. [3] (Man'ko et al. claim that the uncertainty relations are 
necessary to ensure positivity: they are again right, of course, but they do not 
prove this fundamental fact!). Secondly (Section[3]) we propose a new criterion 
for deciding when a phase-space function which is dominated at infinity by a 
phase-space Gaussian is the Wigner distribution of a mixed (non necessarily 
Gaussian!) state. Our approach is based on Hardy's uncertainty principle [5] 
for a function and its Fourier transform. Hardy's theorem goes back to 1933: it 
is unfortunate that its usefulness in quantum mechanics has apparently not been 
noticed before! We conclude our discussion in SectionH]) by linking our results to 
the notions of "quantum blob" and "admissible ellipsoid" introduced by the first 
author in [5J [7J [5] , and which provides a canonically invariant notion of phase- 
space coarse-graining, which seems promising in various aspects of phase-space 
quantization. 

Notation. We work in N degrees of freedom; the coordinates of position vector 
x are xi, ...,xn and those of the momentum vector p are pi, ...,pn- Writing x 

and p as column vectors we set z — (^j ■ We denote by o~{z, z') the symplcctic 

product: by definition a(z, z') = (z') T Jz = p ■ x' — p' ■ x where where J = 



symplectic if and only if S T JS = SJS T = J. 

2 Canonical formulation of the uncertainty prin- 
ciple and positivity 

Let p be a self-adjoint trace-class operator on L 2 (M. N ); we have 



where the kernel K satisfies K (x, x') — K(x', x) and is square integrable on 
x R N . The Wigner distribution of p is the real function 




is the standard symplectic matrix. A 2N x 2N real matrix S is 




(1) 



e * 



*™K(x+\y,x- \y)d N y. 



(2) 
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(see Littlejohn [9] for details; our choice of normalization is consistent with that 
adopted in Weyl calculus, but it differs from that in [3] even in the case K = 1). 
We have 

Tr(p) = J K(x, x)d N x = J W{z)d 2N z. (3) 

Assume now that in addition Tr(p) = 1; if p is positive then it is called a density 
matrix. If A and B are two essentially self-adjoint operators defined on some 
common dense subset of L 2 (WL N ) the covariance of the pair (A, B) with respect 
to p is by definition 

A(A,B)p=l(AB + BA)p-(A)~(B)~ (4) 

where {A)~ = Tr(Ap), and so on. Choosing in particular for A the position 
operator Xj = Xj and for B the momentum operator Pj = —ihd/dxj the 
covariance matrix is the symmetric 2N x 2N matrix 

^»-[A(P,X)? A(PP) P ) W 

where A(X,X) p = (A(Xj, X fc )?)i<i,k<iv. A(X,P)p = (A(X V P k )p)i<j, k <N 
and so on (we assume that all second moments exist; this condition is satis- 
fied for instance if (1 + |z| 2 )p is absolutely integrable). The covariance matrix 
is a fundamental object in both classical and quantum statistical mechanics 
because it incorporates the correlations between the considered variables. The 
Robertson-Schrodinger uncertainty principle says that 

(AA-)|(AP,)| > A(Xj,Pj)p+ \T? ,l<j<N (6) 
(AX,)%AP k % > A{X h P k % for j ? k. (7) 

where A(X,P) is the covariance of the pair (Xj,Pj) (see the original articles 
[10l[TT] and the historical discussion by Trifonov and Donev [12]; for a "modern" 
proof in the general case of non-commuting observables the reader could consult 
Messiah's classical treatise [IB"]). 

One has the following fundamental result well-known in quantum optics, and 
used in the study of entanglement and separability (see for instance [T4l I15j ): 

(I) If the self-adjoint trace-class operator p with Tr(p) = 1 is 
positive, that is if it is a density matrix, then the Hermitian matrix 
+ i^-J is positive semi-definite: 



Sp + yJ>0. I M 



(That Sp + ^ J is Hermitian results from the symmetry of and the fact that 
J T = — J). The relation between property (I) and the Robertson-Schrodinger 
uncertainty principle is the following: 

(II) Condition {Sp is equivalent to the Robertson-Schrodinger in- 
equalities @)-(7p. 



3 



That the formulation |5J of the uncertainty principle is invariant under linear 
canonical transformations follows at once from the fact that S is a symplectic 
matrix; this makes the superiority of this formulation on the usual one: it 
replaces the quite complicated and tedious verification of the inequalities 
([7]) by the calculation of a set of eigenvalues. To see why, let us begin by 
giving two definitions. Let M be any real positive-definite 2N x 2N matrix. 
Since JM is equivalent to the antisymmetric matrix M 1 / 2 JM 1 / 2 its eigenvalues 
are of the type ±ipj (j = 1,...,N) with p,j > 0. Ordering the fij so that 
A*i ^ ("2 > • • • > U>n we call the sequence (p,i, ...,/ijv) the symplectic spectrum 
of M; the number p = pi is called the Williamson invariant of M. Williamson 
[T7] has proved that that there exists a 2N x 2N symplectic matrix such that 
M = S T DS with 

D= \0 a) ' A = dia s(^i' ->Mjv) (9) 

( "Williamson diagonal form" ) . Now, one proves [21 [T5] (see [5] for a detailed 
exposition), that condition ([8]) (and hence the Robertson-Schrodinger inequal- 
ities) is equivalent to: 

(III) The Williamson invariant p of the covariance matrix Ep 
satisfies p > ^h. 

We emphasize that the equivalent conditions (I)— (III) are not sufficient to 
ensure positivity; an illustration is the example of Narcowich and O'Connell |4j 
mentioned in the Introduction. It goes as follows: let the function W(x,p) be 
determined by its Fourier transform via 

f e i{ > xx ' + ^W{x',p')dp'dx' = (1 - \ax 2 - ±/?p 2 ) e - (Q2a;4+/3V) (10) 

(with a, j3 > 0). It is easily checked that W is real and that the corresponding 
operator p satisfies Tr(jo) = 1. Narcowich and O'Connell then show that the 
uncertainty principle is satisfied as soon as a and /3 are chosen such that a (5 > 
h 2 /A. However, even with that choice, the operator p is never non-negative 
because the average of p 4 is in all cases given by 

J p 4 W(x,p)dxdp = -24a 2 < 0; (11) 

p can thus not be the density matrix of any quantum state. 

The considerations above explain the difficulties with positivity questions 
occurring when one rescales Wigner distributions as Man'ko et al do in [3]. Let 
in fact z a denote any of the components of the vector z — (xi, ...,xn',Pi, ■■■,Pn) 
and let Z a , Zp be the operators corresponding to z a , zp. Defining the rescaled 
Wigner distribution W x by 

W x (z) = \- 2N W{\z) , A > (12) 
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we have 

' W x (z)d 2N z = I W(z)d 2N z = 1 (13) 



and a straightforward calculation shows that 

1 



A(Z Q , Z p )p = ^A(Z Q , Zp) f (14) 
where p A is the operator corresponding to W , hence 

Let n and /x A be the Williamson invariants of the matrices Sp and Yj-^x , re- 
spectively. In view of condition (III) in last section, we must have fi > If 
now p x is also to be a density operator we must have p x > ih as well, that is, 
equivalently p > -^A~ 2 h. This requires that A < 1. 

In the Introduction section of this Letter we referred to necessary and suf- 
ficient conditions for a self-adjoint trace-class operator to be positive. In fact, 
Kastler [18] and Loupias and Miracle-Sole [T9J, [20] have shown that the op- 
erator p is positive (and hence a density operator) if and only the following 
so-called "KLM conditions" hold: for every integer m = 1,2,... the complex 
matrix F = (F jk (zj, Zk))i<j,k<m with 

F jk ( Zj ,z k ) = e^ a{z ^ Zk) F„W{z 3 - z k ) (16) 
is positive semi-definite; here 

T a W p (z) = J e l ^ z ' z "^W(z')d N z' (17) 

is the symplectic Fourier transform of the Wigner distribution p. 

In [IB], Lemma 2.1, Narcowich shows that the KLM conditions imply that 
Sp + > 0. In view of (I), (II) above we thus have the following necessary 
condition for p to be positive: 

(IV) Assume that the operator p is positive; then its covariance 
matrix must satisfy condition f3|] or, equivalently, the Schrddinger- 
Robertson inequalities 

This statement thus fully justifies and completes the statement of Man'ko et 
al. [5] that "...the uncertainty principle does not determine the quantum state". 

3 Gaussian estimates and Hardy's theorem 

We ask the following question: 

11 Under which conditions on M can a function W such that 
W(z) < Ce~n Mz ' z be the Wigner distribution of some mixed sate?" . 



5 



The answer to that question is given by following theorem, the proof of 
which relies on the following old result due to Hardy [5]: assume that the 
square-integrable function ip and its Fourier transform 



are such that \ip(x)\ < Ce-^\ x \ 2 and \Fijj(p)\ < Ce-^ 2 (a, b > 0). Then 
we must have ab < 1, and if ab = 1 then tp(x) = Ae^^ x]> for some complex 
constant A. If ab > 1 then ip = 0. 

Theorem 1 Let p be a density operator and assume that its Wigner distribution 
W satisfies an estimate 

W{z) < Ce~i Mz - z (18) 

where C > is some constant and M — M T > 0. Then the Williamson 
invariant [i of M must satisfy p,\ < 1. Equivalently: the matrix Sp = -|M _1 
must satisfy the uncertainty principle: £p + J > 0. 

To prove this we will use the fact that there exists an orthonormal system 
of vectors (ipj)j in L 2 (R N ) such that 

W(«) = ^a J -WVi(«) (19) 

3 

with ay = 1, ctj > (see e.g. [5] and the references therein). Let S be a 
symplectic matrix such that M = S T DS with D as in ©; then the inequality 
(fT5|) is equivalent to 



WpOS -1 *) < Ce-i Ax - x e~i Ap - p . (20) 
Integrating successively with respect to the variables pj and Xj we get 

W{S~ 1 z)dp < C A e-^ Ax - x , [ W{S~ 1 z)d N x < C A e-i Ap - p (21) 



with C\ = C f e h Ax ' x d N x. We next observe that 

W{S- l z) =Y,^W^ 3 {S- 1 z) =J2^W(S^)(z) (22) 
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where S is any of the two metaplectic operators corresponding to S (see for 
instance 8, 9 and the references therein). Taking into account the formulae 

W(SiPj)(z)d N p = \S^(x)\ 2 , J W(SiP J )(z)d N x = \F(S^)( P )\ 2 (23) 

it follows that we have 

]>>|^(x)| 2 <C<Ae-* A - , ^a,|F(^)b)| 2 <C A e-* A ^ (24) 
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and hence, in particular, 

\Sip j {x)\<C jA e-^ Kx - x , \F(S^)( P )\ < C^ac-^p-p (25) 

with Cj,A = \/C\/aj. Since A = diag(/ii, /ijv) with /ii > • • • > hn it follows 
that 

\S^(x)\<C^ A e-^ Ax - x , \F(S^)(p)\ < C jA e-^ A P-P. (26) 
We claim that these inequalities can only hold if [i\ < 1. Set 

<Pj{xi) = sy>j(xi,o, ...,o). 

By the first inequality (f2"6")l we have 

|^(^i)|<Q,A e -^ ? . (27) 

Denoting by Pi the Fourier transform in the variable x\ a straightforward cal- 
culation shows that 

J F(Si>j)(p)dp2 ■■■d PN = (2 7 r/l) (Ar - 1)/2 F 1 J (p 1 ). (28) 
and hence, in view of the second inequality (|26p . 

\FiMPi)\ < {£k) {N ~ 1)/2 C iA J e-^f- "^dpa ■ • • dp N (29) 

that is 

|.Fi&(pi)(pi)| < A e~^ (30) 
for a new constant C^a- Applying Hardy's theorem to (|27p and (|30p we must 
have /if < 1, which proves our claim. 

The result above shows the reason for which the rescaling can be used to pro- 
duce negative operators from a positive one: assume that W p (z) < Ce~^ Mz ' z ; 
then 

\- 2N W p (\z) < \- 2N C x e-i M ^ z (31) 

with C\ = \~ 2N C and M\ = \ 2 M. The Williamson invariant of M\ is A 2 /ii 
and the condition A 2 /ii < 1 will be violated as soon as we choose A > l/y^/Jj". 
(In [3] Man'ko e£ work in units in which ft = 1; it is therefore not immediately 
obvious that the procedure they implement to construct non-positive operators 
by rescaling coordinates in a non-symplectic way is tantamount to increasing 
the value of Planck's constant so that the uncertainty principle is violated.) In 
fact, one immediately understands why the rescaling procedure of Man'ko et al. 
works: it consists (in the example they consider) in replacing M by a matrix 
that is "too small" . In addition, as a by-product of our result we recover the 
property that the support of a Wigner distribution can never be bounded in 
phase space. 

Remark 2 Theorem Q] also allows us to recover in a simple way a result of 
Folland and Sitaram J25^ : a Wigner distribution can never be compactly sup- 
ported. Suppose indeed that there exists some R > such that W p (z) = for 
\z\ > R. For any given fi we can always choose C > large enough so that 
W p (z) < Ce - £' z ' for all z. Choosing fi > 1 this contradicts the theorem.. 
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4 Relation with Quantum Blobs 



In recent previous work [6l [7] one of us has introduced the notion of "quantum 
blob" and of "admissible ellipsoid" in connection with the study of a coordinate- 
free formulation of the uncertainty principle. A quantum blob is the image of 
a phase-space ball with radius \fh by a (linear or afhne) symplectic translation. 
An admissible ellipsoid is a phase-space ellipsoid containing a quantum blob. 
Characteristic properties are: 

• The section of a quantum blob by any plane through its center which is 
parallel to a plane of conjugate coordinates Xj,pj has area ^h; 

• A phase-space ellipsoid is admissible if and only if its section by any plane 
through its center which is parallel to a plane of conjugate coordinates 
Xj,pj has area at least ih. 

Moreover: 

(V) An ellipsoid Bm : Mz ■ z < h is admissible if and only 
if c(Bm) > \h, c any symplectic capacity [5T] on M. 2N , and this 
condition is equivalent to £ + % -^J > with £ = |M _1 . 

and 

(VI) The symplectic capacity of Bm '■ Mz ■ z < ft is c(Bm) 

■ , where (pi, ...,fijsr) is the symplectic spectrum of Bm (see [51 

We can thus re-express Theorem [T] in the following coordinate-free form: 

Theorem 3 Assume that the Wigner distribution of a density operator p is 
such that W(z) < Ce _ « Mz ' z . Then c(Bm) > \h where Bm is the ellipsoid 
Mz-z<h 

It can be interpreted in a very visual way as follows: assume that we have 
coarse-grained phase space by quantum blobs S(B(s/K)). Then the Wigner 
ellipsoid of a density operator cannot be arbitrarily small, but must contain 
such a quantum blob. Equivalently: the Wigner ellipsoid must be defined on 
the "quantum phase-space" consisting of all parts of R 2N containing a quantum 
blob. 

5 Conclusion and Comments 

A "simple" characterization of positivity for trace-class operators is still to be 
found. We have given one such characterization for a particular class of pu- 
tative Wigner distributions (those dominated by a phase-space Gaussian). In 
the general case possibly the phase-space techniques and concepts (symplectic 
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capacities) developed in [51 [7] could provide further insight about what such a 
condition could be (cf. Theorem [3] above). The methods proposed in Bohm and 
Hiley [55] could perhaps shed some light on the question; also see Bracken and 
Wood [23] who introduce the interesting notion of "Groenewold operator" to 
study positivity (but from a slightly different point of view) . 

We finally remark that in the discussed paper [3] Man'ko et al. use the 
notion of quantum fidelity (which, besides, plays an important role in the study 
of Loschmidt echo) to prove that the Wigner function of the first excited state 
of the oscillator does not lead to a positive operator when rescaled; perhaps 
their idea could be exploited in a more general context to shed some light on 
the difficult question of positivity? We will come back to these fundamental 
questions in a forthcoming paper. 
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